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^^v . The fermionic Ising spin glass is analyzed for arbitrary filling and for all temperatures. A selforga- 

O^ ' nized 3-band structure of the model is obtained in the magnetically ordered phase. Deviation from 

^\ \ half filling generates a central nonmagnetic band, which becomes sharply separated at T = by 

T— ( ■ (pseudo)gaps from upper and lower magnetic bands. Replica symmetry breaking effects are derived 

J.J ' for several observables and correlations. They determine the shape of the 3-band DoS, and, for 
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^ ; I. INTRODUCTION 

J^ , Magnetic correlations with frustration in insulating fermionic systems are a key topic of modern condensed matter 

• ■ theory. In disordered systems with frustrated random spin-spin interactions the description is found to be highly non 

j^ ] perturbative, comparable for example with problems encountered in disorder-free Hubbard models. Limits, where 

Cj ■ mean field theories become exact and probably include essential physics, exist, but their exact solutions are hard to 

I ' obtain. These complications are fortunately also linked to the richness of physical phenomena, already incorporated in 

■ O mean field equations. The important experimental tool of doping often leads to quantum phase transitions of various 

C kinds and thus further enhances the structure of phase diagrams. It is the purpose of the present paper to contribute 

to the understanding of the low temperature glassy phase, which reacts strongly to doping before its random magnetic 

order is destroyed in favor of a paramagnetic phase. 

Systems with magnetic interaction often realize saturated magnetic order in the ground state at zero temperature, even 

~~^ ' if the interaction is frustrated and random such that magnetic order consists of randomly oriented frozen-in magnetic 

J^ moments. It is well known that spJB. glass order requires a mean field theory with more than one order parameter or 

even an order parameter functioiJja q(x) with Parisi parameter x, in terms of which complete random order can be 

I ' identified: the support of the subset of order parameters, which do not reach their maximum allowed value at T = 0, 

f^ , shrinks to zero and the Edwards Anderson order parameter reaches its maximum. In case of interacting fermionic 

T— I ■ spins, the magnetic saturation can be increasingly depressed by doping until a breakdown of magnetic order occurs at 

On sufficiently large deviation from half filling. The chemical potential /i as a competing field, which breaks particle-hole 

^^ ' symmetry, and the fernrion concentration i'{fi) are hence important parameters. A famous example for doping effects 

on magnetism is the breakdown of antiferromagnetic order in Hubbard type modelgO. The doped fermionic Ising spin 

glass(/S'G'/), the properties of which are at the center of the present paper, contains a couple of close relationships 

with the Hubbard model. Despite very different techniques employed in the two cases, e.g. the replica formalism for 

the disordered model does not show up in the clean Hubbard model, comparable phenomena in the band structures 

r^ ' appear and will be discussed. l-l 

Q In preceding papers, the density of states of the ISGf was analyzed in detail for half-filling3. A pseudogap was 

O obtained in a solution with infinitely many steps of replica symmetry breaking (RSB). Hence the hierarchy of in- 

^ I finitely many steps of breaking a discrete symmetry resembled the effect of breaking a continuous symmetry, which 

commonly leads to soft modes. In the present case, soft single particle excitation energies were created, irrespective 

of two particle soft modes, for which the single particle DoS appears only as a weight. 

^ Replica symmetry breaking introduces a sort of statistical fluctuation effects in fermionic spin glasses, which determine 

■ - ' the band shape. In particular the softening of gap energies occurs like a onedimensional quantum critical phenomenon. 

This single dimension can be viewed as the replica dimension, while a similar role is played by the time in the dynamic 

mean field theory of the infinite-dimensional Hubbard model. It was however also found for the ISGf , that the finer 

structures induced by symmetry breaking in replica space were directly felt in the quantum dynamic behaviour of the 

subclass of fermionic correlations. 

The description of doping and of general arbitrary filling, which affects size, position, and splitting of gaps, has been 
the goal of the present work. The results of this paper are obtained by entangled (rather than parallel) numerical 
and analytical analysis. This 'numerico-analytical' study is based on a replicated field theoretical treatment of the 
random Ising interaction problem. 



Doping and arbitrary filling enforces occupation of nonmagnetic states at T = 0, even where magnetic order is strongly 
preferred. The present paper involves and revisits earlier results for the tricritical phase diagrarriQ away from half- 
filling. 

The low-temperature limit of the non-half filled model and the way the phase transition into the paramagnetic phase 
takes place, remained an open problem up to now. This was related to several problems: first, information about 
the full replica-broken solution (removing a negative Almeida-Thouless (AT) eigenvalue) was only available at half- 
filling, and secondly another AT eigenvalue turns complex due to the replica limit and the question of stability is 
raised again. 

In this article, however, we suppose that standard replica symmetry breaking is sufficient to describe the ordered 
phase, (see a more detailed discussion of the stability properties in Refs. ||,^. 

After the presentation of the model and a short outline of the calculation, we present first the results obtained within 
the replica-symmetric approximation. The three band structure in the low-temperature regime is already obtained 
in this lowest order approximation within a certain range of chemical potentials. The formulas should be easily 
understandable and helpful for any reader, who wants to limit his interest to the selforganized appearence of the non- 
magnetic third band; considerably increased efforts by one step RSB lead to a refined picture of this 3-band structure. 
We discuss the way these three bands emerge below the freezing temperature and how they can be understood in 
terms of magnetic and nonmagnetic contributions. It is well known that replica symmetry breaking occurs and, 
reemphasizing that it has a profound effect on the low temperature properties of the system, we show how it is man- 
ifested together with broken particle hole symmetry. For this purpose we present in detail calculations and results in 
one-step replica symmetry breaking (IRSB). Fortunately, this solution already allows a sound estimation of the exact 
one. Some key features suggested by this one-step RSB solution are also confirmed by exact results for infinite RSB: 
among those, the replacement of the hard gap by a pseudo gap and the onset of the three-band splitting already at 
arbitrary small deviations from half-filling are most important. 



II. THE DOPED FERMIONIC ISING SPIN GLASS 

Since in several preceding publications the model has been explained in many respects, we wish to limit the present 
discussion to the specific aim of this paper, hence the way doping interferes in the interplay of charge- and spin- 
correlations. 



The grandcanonical Hamilton operator 



^ = Z! ^y ^'^l ~ /^ Z! '^''^ (^) 



i<] 



describes instantaneous interactions of fermionic spins. Stripping off the factor h/2 the spin operators a^ are given 
by the fermion particle number operators by cr^ = n| — n^ . The distribution of random interaction couplings 



generates magnetic correlations independent of time and of infinite range in real space. 

III. MAGNETIC AND NONMAGNETIC BANDS 

A. Replica Symmetric Results 

Despite its instability against replica symmetry breaking (RSB) the replica symmetric solution is nontrivial and 
must be understood in detail, since it forms the basis for the improved solution presented below. Moreover, in a 
relatively simple description, this approximation contains features of the selforganized transition from single-band 
structure above freezing temperature to either two-bands for half filling or three-band structure, which occur below 
freezing of the magnetic moments. The quality of this approximation decreases in the low temperature regime, as 
the next section will show. The role of symmetry breaking amounts to a softening of the gaps and to a remcojal of 
sharp drop-offs of the density of states at low temperatures, as demonstrated earlier for the case of half-fillingP. At 
half-filling a pseudogap results between upper and lower magnetic band. A special line in the phase diagram windiiic 
through the ordered phase, connecting tricritical point(s) and the zero temperature gap edges of the magnetic bandsH, 



wraps a regime where the spectral weight of the central nonmagnetic band does not fully move into the magnetic 
bands as T ^- 0. 

In this section we first present selfconsistent numerical solutions for the temperature range from the filling-dependent 
freezing temperature down to very low T of order 10^^ J. The numerical evaluation of the density of states employs 
numerical solutions of the selfconsistent equations for spin glass order parameter q and linear susceptibility % as a 
function of temperature and parametrized by the chemical potential /i. The fermion concentration is also calculated 
as a function of /i. This analysis is supplemented by an exact T — Q calculation. 

1. The free energy and resulting selfconsistent equations at finite temperatures 

The following results are derived from the effective Lagrangian and from a generating functional for general fermionic 
correlation functions of the model. Analytical solutions are found in the sense that the number of nested integrations 
becomes minimized before, in a final step, the observable are evaluated numerically. In order to facilitate the T — )■ 
limit, it is useful to change variables by introducing the linear susceptibility x = PJifl ~ <?)■ Since q — q decays 
linearly with T in the T — > limit, the finite linear susceptibility is a helpful quantity which reduces the degree of 
/?-divergences (which need to be compensated). For this reason / is expressed rather in terms of q and x than in 
terms of q and q. 
The free energy (density) / is obtained in the replica symmetric approximation as 

f^^Jx{jj + 2q-2)-Tln2-ti-Tj InC with (3) 

C == cosh(/3jygz) + co8h{Pti)e-i^-'^ (4) 

The Gaussian integral, which appears frequently throughout the paper, has been reduced to the short form defined 

by 

0(z)EE— =/ dze-^V(^) (5) 



Extremalization of expression (0) with respect to q, Xi ^-^d /i yields the coupled selfconsistent equations 

= dJ^dJ^iy + d^f (6) 

which we solved numerically. The results are finally used in the band structure calculation. 

2. Energy and selfconsistent equations in the T — limit 

The zero temperature limit deserves separate attention for several reasons. As can be seen from the thermal 
free energy, cancellation of divergences in the /? — > oo limit trouble the numerical work at low temperatures but 
nondivergent analytical results at T = provide a control. 

For T = 0, an analytical approximation was also obtained, applying an expansion in powers of /i — x/2. Good 
agreement was found almost up to the magnetic breakdown. ■— ■ 

As usual the equations are simplified in the T = limit and were obtained by a variantcJ of the Sommerfeld method. 
For the density of states the steepest descent method is used in addition. 
We obtain for /i < ^Jx'- 

£; ^ /(T = 0) = ij x(g - 1) - A^ - ^ J^/9 (7) 

while for /i > 5 JX; the T = energy becomes 

!,,(„ II .. (.. (,,1 -I- I 1 11 I.I „<1V 



E = -Jxi, -I) -,-{,- -^JX) erf (^(5 - |)j - J^ ^^^^.-^^ (8) 



Extremalization of these energies yields T = self-consistency equations which couple the magnetic correlations q 

and X with the filling factor as a charge average. 

For J > f we derive the following relations between zero temperature parameters 



"^' + "'145(7 -I')- 



X = 



e 29 




^(f-f)' 



(9) 
(10) 



One may derive the T = solutions as a function of either v or /i. The T = relations q = q and q — 2 — ly hold also 
for < 4 < ^, whence in this interval one simply obtains 



(11) 



3. Density of states 



The derivation of the Green's function from the generating functional of spin glasses was discussed beforeQ. The 
present work makes use of the same formalism, but evaluates it for all fillings. The (numerical) solutions for q{^i,T) 
and xiPi ^) ^re employed in the calculation of the electronic density of states throughout the whole spin glass phase. 
The following set of Figures |l] shows that a central band emerges for high enough chemical potentials. In the ORSB 
approximation the band gap discussed previously for half fiUingQ is visible up to fi = ^Eg — --h= and the DoS looks 
like one half of a bath tub. When the chemical potential is moved towards the band gap value (we show for example 
the DoS dA. fj, — 0.39), a tiny central band shows up at low temperatures, but looses its weight again completely in the 
T -^ limit in favor of the magnetic side bands. Roughly speaking, a line given by d^f/dq^=0, which, starting at the 
tricritical point of the phase diagram, bends into the gap edgeB at T = 0, wraps this small precursor of the central 
band. For chemical potentials exceeding the gap edge value -h=, the central band manifests itself already at higher 
T and survives at T = 0, becoming there a pure Gaussian function oi E = e + fi with finite height and symmetric 
cutoffs. In the ORSB approximation both the central band and the magnetic bands are sharply cut off and separated 
by gaps of identical width x- As replica symmetry breaking will show, the gap size is given by the nonequilibrium 
susceptibility x, which agrees with x only in this lowest order approximation - we continue to discuss the gap in this 
chapter in terms of x- This quantity, which separates central from upper band and central from lower magnetic band 
begins to vary with /i for |/i| > 
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FIG. 1. Density of states p{E — e + iJ.,T) as a function of temperature below freezing, as calculated in the replica symmetric 
approximation at specific chemical potentials /i (in units of J). The bold front line represents the independent T = result 
obtained from Equations (7-12). 

Figure m combines the numerical finite temperature calculations of p{e,T\n) with the exactly calculated function 
p(e, T = 0|/i) given below, making use of the numerical solutions for order parameter and susceptibility. As emphasized 
by the thick lines at T = 0, the T — density of states can be decomposed into isolated contributions of the three 
bands. Taking advantage of the symmetry one may use the energy variable E = e + j.i 



piT^O,E)=p^iO,E) + p,{0,E)+p+iO,E) 



(12) 



where the central charge band is given by 



PciO,E) = 



^2^9(0,^) 



.-E^/i2,io,,y) e(A.-x(0,A^)/2-|i?|) 



Upper and lower magnetic band contribute 



P±{E)^ 



1 



^2^9(0,^) 



(e-(|s|-x(o,M))V(2«(o,,0)e(|ij| _ (^ + ;^(0,/,)/2) 



(13) 



(14) 



Let us first discuss the numerical results in the RS-approximation. Figure ^ shows a collection at typical values for 
the chemical potential: 

i) Within the hard gap regime < |/i| < l/V^n, where v{T = 0) = 1, a pronounced central band is absent; only 
in a small range of low but finite temperatures a tiny midgap peak is observed. We find that its existence is clearly 
linked to the characteristic line mentioned above. This line separates the domain of the phase diagram, where the free 
energy is minimized as a function of x, from the one where it is maximized. This latter property is unrelated to the 
wellknown maximization of / by the SG order parameter q; in turn the presence of q is needed to render a solution 
with d^f/dx^ < stable. For the present case we find a complex Almeida Thouless eigenvalue (for replica-diagonal 
perturbation), which does at least not exclude stability apart from the Parisi RSB. 

ii) For chemical potentials large enough to sustain fillings different from one electron per site, i.e. |^| > — i=, three 
bands develop as the temperature falls below Tf and become separated as T ^ 0. The Fermi level lies between the 
central charge band and the upper magnetic band. The area under the central peak belongs to the deviation from 
half-filling, described by 1 1^ — 1 1 . The low T results (see Figure ||) confirm the numerical observation that the central 

band width is given by E^^^' = 2/i — x at T = 0, that both left and right gap widths obey Eg ' = X; while at half 

filling the relation reads Eg ' — 2%. 

One may compare this approximate solution with an iterated perturbation solution of the half-filled Hubbard model 

at zero temperatureEl: there a hopping generated band shows up within the Hubbard gap insulating phase. 




FIG. 2. The density of states shown for temperatures T — .3, .1, .01 at chemical potentials fi — .1, .39, .41, .5, .6, .7, .8, .9 




FIG. 3. DoS at zero temperature as a function of /x. Above fi = I/vStt, the central charge band appears. In this case, the 
Fermi level lies between the central and the upper magnetic band. 



B. Improved solutions with broken replica symmetry and strong low— temperature effects 



One step replica symmetry breaking yields a large step towards the exact solution: it allows to guess properties of 
the exact one quite frequently. As the basic starting formula for all thermodynamic properties we rewrite the free 
energy in IRSB in terms of parameters, which also allow to obtain finite zero-temperature limits. It is therefore 
convenient to use the nonequilibrium susceptibility, linear susceptibility, and Parisi parameter denoted respectively 
by X = PJ{q ~ 9i)) X = X + '^('Zi ~ 12), and a = j3Jm (we have chosen x, X) and a as dimensionless) . 



Now the free energy density reads 



1 



^ ^i-^(^(#7 + 291 - 2) + ci{ql ~ 92)) - Tin 2 - /i - - 
C = cosh(/3 J(V9i - q2Zi + ^/g2Z2)) + cosh(/3/i)e~~^ 



In 



C^ with 



(15) 
(16) 



The condition of stationary free energy provides saddle-point equations for the four parameters x, qi , g2 , and a, which 
need to be determined as functions of fi and T. The derivation of the thermal selfconsistent equations is lengthy and 
hence omitted. Instead we present in Figure ^ the solutions as a function of /i for three selected values of temperature. 
For the use in the density of states below, we have determined all necessary parameters almost continuously on a grid 
of (A^/J, AT/J) = (10-2,10-2). 
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FIG. 4. Order parameters, nonequilibrium susceptibility, and Parisi parameter, shown in IRSB, as a function of // at selected 
temperatures T — 0.05 (solid lines), T = 0.15 (dashed lines), and T = 0.3 (dotted lines). 



The above formulation of the free energy allows to obtain the zero temperature limit in terms of finite quantities 

(17) 



/ = 7;Jx{qi - 1) + -rJaiqi ~ql)~ fj- — 
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The zero temperature limit of the internal integral is calculated as 
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(18) 



(19) 



These preceding T = 0-limits are the basic ingredients which we use analytically to get the selfconsistent equations 
at zero temperature in terms of finite parameters. Results of our final numerical evaluation are collected in Figure 0. 
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FIG. 5. Order parameters in IRSB at zero temperature 



All of these parameters show a remarkable variation before the first order transition regime to the paramagnetic 
state is approached at ^j,ci{T = 0) « 0.881. Even gi, which agrees with the spin autocorrelation function q at T = 0, 
contains interesting behaviour. This will be extracted below in terms of the change of the fermion filling under one- 
step replica symmetry breaking. 

The filling factor is an integrated quantity over the density of states folded with the Fermi distribution. We can also 
exploit all solutions, given so far for the order parameters, in order to determine the fermionic density of states itself. 
This is done first for all relevant chemical potentials and all temperatures O(10~^) < T < Tf{iJ.), resulting in the 
body of the Figures Q The T ~ solutions are then derived as follows. Matching perfectly the results at lowest finite 
T, they arc finally combined with the body of each figure. 



The density of states reads in IRSB for T > and as a function oi E = e + fi 

rG 



T I a gj_ rG f'^ C'^-'^e.-P(-E+(j^q:^z^+.j^^T^^zt)Y /{2Jx) 

p{E) = ^J^(cosh(/3/i) +cosh(/3£:))e~-'^ ' ^^ 
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^\ Jx 

In the T = limit, the saddle point method allows to solve the internal integrals exactly, which results in 
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for ^ < J| and \E\ < Jx 
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for /i > J| and \E\ < ^ - 

Jx 



2 



for /I > J| and fi ~ ^ < \E\ < ^i + ^ 



Jx 






for /i > J^ and /i 



Jx 



<\E\ 



where Eqs. (nsf) and dl9) have to be substituted. 
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FIG. 6. Density of states at finite and at zero temperature in IRSB shown for fi = .05, 
E = e + fi and temperature T. 



and .7 as a function of energy 



Figure ra displays the exact selfconsistent evaluation at one step replica symmetry breaking. The analytical inter- 
mediate solutions described before, allowed to include the T = solution into the Figure. The full picture shows the 
band structure at all temperatures below Tf. The choice of fi values covers a wide range from almost hal filling at 
H = 0.05 to /x = 0.7 through the entire spin glass phase. 

A comparison with the lowest order approximation, shown in Figure |l|, combined with analytical results for ooRSB, 
gives a clear hint for the exact ISGf solution. 

Only the first part at /i = 0.05 < ^Eg ' ss .119 does not contain the central band, since j/(/i = .05, T = 0) = 1 

in IRSB. For E^ ' = .119 < fi < Eg — J ^ the central band is present in IRSB, while it was absent in ORSB 

in this interval. This effect sets in at low temperatures, corresponding to the smaller magnetic energy scales set by 
the new additional order parameter of IRSB. Above temperatures determined roughly by the random field condition 
d^ f /dq"'"''^ — Q {q = q°'°'), which corresponds to a random field crossover line Trf{fJ,), the spectral weight does not 
show a peak at E — 0. It needs larger values of /z to find a large central peak already for intermediate temperatures 
below Tf. On the other hand, the central band is already well developed at /z = 0.5 for low enough temperatures. At 
this value of /i the lowest order approximation leads only to a very small band, since /i — 0.5 exceeds only by little 
the ORSB gap edge. 

As for half filling the ratio between the gap widths and the (finite) DoS at the gap edge is constant. For higher order 
symmetry breaking the gap shrinks and the spectral weight at the edge diminishes correspondingly; thus spectral 
weight is moved into parts of the gap region, as the approximation is improved step by step. 

The crossover from finite low T to the exact T = solutions, shown by fat lines in the 3D-plots, redistributes consid- 
erably the spectral weight. At higher /u, for example /i = 0.7 as shown in Figure 0, the central band shows a maximum 
also as a function of temperature. The DOS height in the center decreases again at lowest temperatures, the band 
becomes broader at T = 0. 

Comparing the T = limit of Figures 1 and 2, which show the ORSB approximation, with that of Figure || for the first 
improved IRSB approximation, one finds that replica symmetry breaking leads to a refinement of the band structure. 



Each further RSB step, until the exact ooRSB sohition is reached, will modify magnetic and nonmagnetic bands 
according to the still smaller scales and magnetic order parameters. However, although one step RSB is not yet exact 
in this model, further refinements are much smaller in size and it is possible to imagine the exact result from our 
Figure H- 

As a function of a continuously varying chemical potential, the IRSB solution for the zero temperature density of 
states is displayed in Figure 0. 




2^0.0 
FIG. 7. DoS at zero temperature for IRSB. The central band appears above /i ~ 0.119. 

The central band maintains the wedge- like shape already observed in the ORSB result. Its thick end however shows 
a new structure near the discontinuous magnetic breakdown, while the other end progressed by a large step towards 
/i = 0. The distance from /i = will be further reduced in kRSB with k > 1. 

1. Infinite breaking of replica symmetry 

The top of the wedge reaches ^ = in ooRSB and at the same time the gap widths approach zero; the DoS becomes 
zero at ±/i but stays finite yet very small near E = zb/u. The derivation is analogous to the one presented in Ref. [g: 
The ratio between gap width DoS value at the gap edge is invariant under RSB and hence they decay together to 
zero in the limit of infinite RSB. Here it is assumed that the nonequilibrium susceptibility, which determines the gap 
width at any order k of the RSB, approaches zeroJor fc -^ cx) as in the half- filled case. The latter result was inferred 
from the work of Thoulcss, Anderson, and Palmcrtil. 

IV. THE CONNECTION BETWEEN FERMION CONCENTRATION AND REPLICA SYMMETRY 

BREAKING 

The filling factor can be described by the summation over the imaginary frequency Green's function Q by 

j. = T^g,(e„)e-"+ (22) 

The Green's function is related to the density of states discussed before by the usual spectral representation 

t/(e;) = J de Z' _ (and by p{e) = — -/m[G^(e)]), which means that the whole fermion propagator changes under 

RSB. Still the summation over all frequencies could either wipe out or maintain this dependence. Indeed what we 

find is a transition between these two alternatives inside the ordered phase. 

We have emphasized the role of replica symmetry breaking for quantum dynamics and low energy excitations. Despite 

the absence of quantum dynamics in charge correlations, and the absence of spin-charge couplings in the Hamiltonian, 

we find that replica symmetry breaking affects spin-correlations and charge-average in a qualitatively different way. 

This is remarkable also because i/ is a global quantity (summed over all frequencies) in contrast to the Green's function 

(which is local in e). 

In one-step breaking we report in this section a crossover line, which separates a regime of almost invisible RSB-effects 

in the filling factor and in q from one with large RSB-effects in these quantities. The magnetic order parameter does 

10 



not show any sign of this crossover. The announcement of these effects evokes on one hand the old, resolved problerntHl 

of absence of RSB below the freezing temperature, and on the other hand thei^henomenon of a Gabay- Toulouse linetJ, 

followed by a crossover to a region with RSB effects in all order parameterstil. 

The first case bears no relation |W|ith our case: the present model is static in charge and spin-correlation and the 

particular role of dynamic effectst3, which occurcd in the transverse field Ising model, do not exist here. 

The crossover line describing the onset of RSB in transversal correlations of a Heisenberg spin glass in a magnetic 

field however, has a vague resemblence, provided we imagine charge degrees of freedom as transversal with respect to 

spin. The chemical potential then roughly corresponds to the magnetic field in the standard case. However a detailed 

mapping between the two models does not seem feasible. 

In Figure || the fermion filling factor is shown. Pairs of {i'orsb{iJ.) ii^irsbifJ-)) are grouped together for ii = .1, .2, ..., .8. 

The detailed plots contain two interesting features: each pair of lines seems to merge asymptotically, but the lines 

still cross each other at a Tq (/i), staying close together for Tq (^) < T < T/(/z). Since the lines cross, having almost 

identical slope, it is difficult to determine with sufficient numerical precision the line Tq (/z) of crossings. We therefore 

chose points T*{fi) where ORSB- and IRSB-lines differ only by 10~^. In between these points and the endpoints at 

the freezing temperature corresponding to the /i-parameter of each curve, the ORSB and IRSB curves cross at least 

once. 

The interesting fact now is that there is a region below the freezing temperature, where replica symmetry breaking 

effects in the fermion filling factor almost vanish. While this occurs in the charge related quantity, the magnetic 

observables such as order parameters still show large RSB effects. One may find it more surprising that RSB effects 

appear at all in the fermion concentration, but it seems still more important that these effects disappear almost (or 

perhaps completely in the exact solution of infinite RSB) within the spin glass phase. 

We do not attempt here to discuss any result for itinerant models, but it is nevertheless clear that the feature described 

for the insulating model can have implications on itinerant systems. 

In addition to the separation of magnetic and nonmagnetic bands described before for the ordered phase, one also 

finds a different exposure of charge- and spin-related quantitites to the effects of replica symmetry breaking. None of 

the quantities is excluded from this, despite the fact that the model does not contain a spin-charge interaction. The 

only coupling is mediated by the chemical potential. 




0.1 0.2 0.3 



0.4 0.5 



FIG. 8. Results for the temperature dependence of the filling factors in ORSB and IRSB are grouped together for .1 < /i < .8 
with A/i = .1. The right endpoints are the freezing temperatures T/ (/i), while at the intermediate points, denoted by T* on the 
pair of lines for /i = .8, the fillings vqrsb and v\rsb differ by only 10"'*. Zeros of vqrsb — virsb exist within T*{fi) < T < Tf{^) 
and at low T for fj, > .7 
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FIG. 9. The filling factor at T = are shown as a function of the chemical potential /i. The ORSB- and IRSB-results of 
this paper are included in the fisare on the left hand side, while on the right the numerical solution of the corresponding 
TAP-equations is taken from RefJij for comparison. There, N denotes the system size. 

Beyond the crossing point of the two approximations at /x* « 0.7 J both hnes stay close together as if no symmetry 
breaking effect would occur at all between /x* and the discontinuous breakdown of magnetic order at /zi « 0.881. On 
the other side, for fi < ji* the RSB effect is large. Moreover, the right hand-^ide of Fig.p shows for comparison the 
numerical result for the corresponding TAP-equation of finite size systemst-3. This generalized TAP-result should 
correspond to the full Parisi solution with ooRSB. At T = the filling i/ starts to differ from one, when the chemical 
potential moves through the gap edge. The gap size depends on the number of RSB steps and decreases to zero in 
the cxjRSB solution. Thus it is clear that left endpoint of the filling curve i/(/i) moves into (i^ = 1, /i = 0). The change 
from the calculated 1-step RSB solution to this exact one is smaller than the one from ORSB to IRSB. The shape 
of the IRSB solution resembles almost perfectly the one found numerically from the generalized TAP-equation. Also 
quantitative agreement is obtained for fi smaller than roughly 0.5. The nonvanishing kRSB-corrections for k — 2 
and higher can be expected to be almost invisibly small. Deviations in the high-/i region might be due to numerical 
problems of the algorithm in na either because of the vicinity of the first order phase transition or due to finite size 
effects. 

In general however the TAP-solution is already in good agreement with the analytical IRSB solution, whose extension 
to the full RSB solution is obvious. The flat increase of z^(/i) from 1 (probably with slope zero) in the generalized 
Parisi solution is consistent with the fact that for rare nonmagnetic regions the central band does not start as a 
5-peak, but as function with finite height and a width increasing smoothly from zero as /x becomes finite. 



V. SUMMARY AND OUTLOOK 

By analyzing the spin glass phase in the (/i, r)-plane we found two magnetic and one central nonmagnetic band, 
which become separated perfectly at T = 0. Their separation by finite gaps in any finite RSB approximation turns 
into a pseudogap separation in the exact infinite-step RSB solution, assuming that the nonequilibrium susceptibility 
vanishes as predicted by the TAP-solution. 

RSB effects did appear in the magnetic order parameters, in the density of states, and in the quantum dynamics 
displayed by the Green's functions, and also in the fermion filling (being at T = simply the integrated DoS). 
The results presented here for the ISGf model trigger speculations not only on other insulating random interaction 
models, like XY- and Heisenberg quantum spin glasses, but also on itinerant extensions. 

For example, the photoconductivity of the itinerant extension of the model is related to the integrated overlap of 
frequency-shifted density of states. The filling factor proved that integration does not wash out RSB-effects, which 
could therefore be expected to exist in the photoconductivity at low temperatures. In itinerant systems the possibility 
of localization in the pseudogap regime of small DoS will be decisive for ac- and dc-conductivity. Again RSB-effects 
will emerge and mark the theory of localization due to a frustrated random magnetic interaction in the spin glass 
phase. It will be highly interesting to compare this with the fully frustrated Hubbard model in infinite dimensionsliS. 
While any finite step RSB does not yet allow to state that the exact line T* (/i) must be different from Tf (/i) (as it is 
in IRSB) and that RSB perhaps exactly vanishes in the region T*(fj,) <T < Tf{ii), this remains a serious option. 
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For itinerant models the question arises whether an RSB transition, disconnected from the magnetic transition, can 
also occur in metallic spin glasses and perhaps will influence transport properties. 
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